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We consider two-dimensional (2D) localized modes in the second-harmonic-generating (x' 2 "*) sys- 
tem with the harmonic-oscillator (HO) trapping potential. In addition to its realization in optics, 
the system describes the mean-field dynamics of mixed atomic-molecular Bose-Einstein conden- 
" sates (BECs). The existence and stability of various modes is determined by their total power, N, 

topological charge, m/2 [m is the intrinsic vorticity of the second-harmonic (SH) field], and 
mismatch, q. The analysis is carried out in a numerical form and, in parallel, by means of the 
variational approximation (VA), which produces results that agree well with numerical findings. 
Below a certain power threshold, N < Nc" l \q), all trapped modes are of the single-color type, 
represented by the SH component only, while the fundamental-frequency (FF) one is absent. In 
1 contrast with the usual situation, where such modes are always unstable, we demonstrate that they 

are stable, for m — 0,1,2 (the mode with m = 1 may be formally considered as a semi-vortex 
with topological charge m/2 = 1/2), at N < Ni m \q), and unstable above this threshold. On the 
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. other hand, TVi" 1 ' (q) = at q > q mxx (in our notation, q max = 1), hence the single-color modes are 

unstable in the latter case. At N = Ni m \ the modes with m = and m = 2 undergo a pitchfork 
0-4 bifurcation, which gives rise to two-color states, which remain completely stable for m — 0. The 

two-color vortices with m = 2 (topological charge 1) have an upper stability border, N = N C 2(q). 
ryj 1 Above the border, they exhibit periodic splittings and recombinations, while keeping their vorticity. 

O ! The semi- vortex does not bifurcate; at N — iVi" 1 : \ it exhibits quasi-chaotic oscillations and a 

ryj ■ rotating "groove" resembling a screw-edge dislocation induced by the semi-integer vorticity. 

x, 

OSIS numbers: 190.6135; 190.3100; 190.4410; 020.475 



I. INTRODUCTION 



It is commonly known that the quadratic, alias x^ 2 \ nonlinearity, plays a fundamental role in nonlinear optics, 
helping, in particular, to create various species of solitons [U-Q, @- The use of the x^ nonlinearity is crucially 
important for the making of 2D and 3D solitons, because, on the contrary to Kerr (x^) nonlinearity, the quadratic 
interaction between the FF and SH fields does not give rise to the collapse Q, which is a severe problem for the 
stability of 2D and 3D solitons in Kerr media Q. Thanks to this circumstance, the first y/ 2 ' solitons were created 
as stable (2+l)-dimensional beams propagating in an SH-generating crystal Further, the absence of the collapse 
instability in the 3D setting suggests a possibility of creating fully localized "light bullets" In the experiment, 

3D solitons have not been observed yet, the best result being a spatiotemporal soliton self-trapped in the longitudinal 
and one transverse directions, due to the interplay of the diffraction, group- velocity dispersion, and nonlinearity, 
while the confinement in the other transverse direction was provided by the waveguiding structure [111 Il2| . 

Another natural possibility in the (2+l)D setting is the creation of vortical solitary beams, with the "hollow" in 
the middle. In these modes, self-trapped SH and FF fields carry intrinsic vorticities m and m/2, respectively. The 
modes are classified as vortex solitons with topological charge m/2; solitons with odd values of m are not possible, as 
the intrinsic vorticity of the FF component, m/2, cannot take half-integer values, although vortices with a half-integer 
optical angular momentum can be created, in the form of mixed screw-edge dislocations, by passing the holding beam 
through a spiral-phase plate displaced off the beam's axis [13j. Unlike their fundamental counterparts with m — 0, 
the vortex solitons in the free space are always unstable against azimuthal perturbations, which split them into sets 
of separating segments. This instability was predicted theoretically [13]- [13 an d demonstrated in the experiment [l8j . 
The same instability was also predicted in the framework of the so-called Type-II (three-wave) y/ 2 ) system, which 
includes two distinct components of the FF field [Til H3| • 

Solitons with embedded vorticity are also known as solutions to the 2D nonlinear Schrodinger (NLS) equation with 
the self- focusing cubic term [2lj ^and they too are subject to the azimuthal instability, which is actually stronger than 
the collapse-induced instability Q • The 2D self- focusing NLS equation models not only the light transmission in bulk 
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media with the Kerr nonlinearity, but also [in the form of the Gross-Pitaevskii (GP) equation] the mean-field dynamics 
of BECs in ultracold gases with attractive inter-atomic interactions, shaped as "pancakes" by the confining potential 
[2^ |. A solution to the instability problem was elaborated in the latter context: both fundamental solitons and 
solitary vortices with topological charge 1 can be stabilized by isotropic HO (harmonic-oscillator) trapping potentials. 
As shown in detail in a number of theoretical works [23j-[29], the HO potential stabilizes the fundamental solitons in 
the entire region of their existence, while vortex solitons are stabilized, in terms of their norm (the counterpart of the 
total power of spatial optical solitons), in ~ 33% of their existence region, and in an adjacent region of width ~ 10% 
vortices exist in the form of periodically splitting and recombing modes, which keep their vorticity [27j . 

The effective 2D trapping potential can be also realized in optical waveguides, in the form of the respective profile 
of the transverse modulation of the local refractive index [j| . This circumstance suggests a natural possibility for the 
stabilization of (2+l)D vortex solitons in the medium by means of the radial HO potential, which is the main 
subject of the present work. A feasible approach to the making of the optical medium combining a nearly-parabolic 
profile of the refractive index and x^ nonlinearity is the use of a 2D photonic crystal, which can be readily designed 
to emulate the required index profile, while the nonlinearity is provided by the poled material (liquid [6( or solid 
0) filling the voids. As shown below, the effective radial potential provides for sufficiently strong localization of the 
trapped modes, therefore the exact parabolic shape of the radial profile is not crucially important. The analysis can 
be readily adjusted to other profiles, if necessary. 

The model, based on the system of coupled equations for the FF and SH fields, is introduced in Section II. It is 
relevant to mention that essentially the same system of GP e quat ions for the atomic and molecular mean-field wave 
functions describes the BEC in the atomic-molecular mixture |30||-[33j|. Accordingly, the predicted mechanism of the 
stabilization of two-component vortex solitons trapped in the HO potential can also be realized in the BEC mixture. 

Solutions for the trapped 2D modes and their stability against perturbations are considered in Sections III and 
IV. First, we address states which, in the unperturbed form, contain only the SH field, while the FF component 
vanishes. Such single-color solutions of the x^ system are known in other contexts, but they are usually subject 
to the parametric instability against small perturbations in the FF component. Our first result is that the trapped 
single-color modes, both fundamental and vortical ones, have a finite stability domain. In particular, the single-color 
modes with m = 1, which (formally) look as semi-vortices with topological charge 1/2, are also found, and they 
are stable too in a finite parameter area. The single-color states become unstable at a particular critical value of 
the total power (alias norm); however, the critical norm vanishes if the mismatch parameter, q, is too large, viz., 
at q > 1 in the notation adopted below). As concerns the fundamental modes (m — 0) and those with topological 
charge 1 (m = 2), at exactly the same critical point they undergo a pitchfork bifurcation, which gives rise to two- 
color states. The stability of the fundamental two-color mode is obvious, while a nontrivial result is finding stability 
borders for the trapped two-color vortex soliton. Above the instability border, it develop periodic oscillations, keeping 
its vorticity and featuring periodic splittings and recoveries. As concerns the single-color semi-vortex, it exhibits a 
different behavior, developing persistent quasi-random oscillations above the stability border, mixing the first and 
zeroth angular harmonics in both the SH and FF components. 

The results concerning the shape and stability of the modes of all the above-mentioned types (single-color ones with 
m = 0,1,2, and two-color complexes with m = and 2) are obtained, in parallel, by means of numerical methods 
and in an analytical form, based on the variational approximation (VA). In almost all the cases, the VA demonstrates 
very good accuracy in comparison with numerical results. 



II. THE MODEL 

The model for the SH generation in 2D is based on the usual scaled equations [H-0 for the FF and SH field 
amplitudes, u and v, 

,du 1 / d 2 d 2 \ , „, , 
l d- Z + 2{d^ + W) U + UV ~ U(X ' V)U = °' 
dv 1 / d 2 d 2 \ 1 



where z is the propagation distance, x and y are the transverse coordinates, the asterisk stands for the complex 
conjugate, q is the real mismatch coefficient, and, as said above, the axisymmetric modulation of the refractive index 
is modeled by the isotropic HO potential, U(x,y) = (fi 2 /2) (x 2 +y 2 ). By means of an obvious rescaling, we fix 
f2 = 1/2, while q remains a free parameter. 

The relation between the potential terms in the equation for the FF and SH fields implies that the same refractive 
index acts on both fields, which is a realistic assumption for materials of which the above-mentioned pair of the 
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photonic crystal and voids-filling stuff can be fabricated. If the weak index dispersion is taken into regards, it will 
produce only a small perturbation in the system. 

The 2D system of scaled GP equations for the atomic-molecular mixture corresponds to Eqs. fl) with z replaced 
by time t, and q accounting for a difference of the chemical potential between the atomic and molecular components, 
u and v. It is relevant to mention that a similar model with a spatially periodic (lattice) potential U (x,y) was 
considered in Ref. 34], where it was demonstrated that the lattice can readily stabilize two-color solitary vortices, 
although with an anisotropic shape. 

Equations (Q]) are derived from the corresponding action, S — j Ldz, with Lagrangian 



L 



iu z ii* + 2iv z v* 



(KI 2 + KI 2 



-U(r) (\u\ 2 + 4M 2 ) - q\v\ 2 + i (uV + u 2 *v) | dxdy. 



(2) 



Stationary modes can be characterized by their total power (norm), N = J J (\u\ 2 + 4|t>| 2 ) dxdy. 



III. STABILITY OF SINGLE-COLOR BEAMS 

A. The beams with topological charge and 1 (m = and m = 2) 

For the single-color modes with u — 0, the SH field satisfies the 2D linear Schrodinger equation with the isotropic 
HO potential, 

dv 1 / d 2 d 2 \ 

*% + - 2 { 1 ^ + w )v-<iv-mx,y)v = o. (3) 

Stationary solutions to Eq. are commonly known from quantum mechanics. In polar coordinates (r,9), they are 

v = v m0 exp (i (-fj,z + mff}) r m exp (-r 2 /2) , (4) 

with arbitrary amplitude v m o (it is defined to be real), integer orbital quantum number m (as said above, it corresponds 
to the beam's topological charge m/2), and eigenvalue /j, = (m + 1 + q) /2 (recall £1 — 1/2 is fixed). The total power 
of this solution is 

N m = ^m\v 2 mQ . (5) 

Solutions (|4]) are obviously stable within the framework of linear equation ([3]), the issue being to find a threshold, 

N = Nc (q), at which the parametric instability against infinitesimal perturbations in the FF field sets in, due to 
the nonlinearity in Eq. (p} for the FF field. As shown in Fig. HJa), the threshold was identified from systematic 
simulations of the perturbed evolution of the single-color beams within the framework of the full system of Eqs. (JlJ, 
for m = and 2. 

The stability was also investigated in an analytical form by means of the VA. For m = 0, the ansatz for the 
perturbed solution is taken as 

u = uq(z) exp (-ar 2 ) , v — vq(z) exp (— r 2 /2) , (6) 

where amplitudes uq(z) and vq(z) are treated as variational parameters, while a is a free constant, which is used 
below as a variational parameter too, but in a different sense. The substitution of ansatz ([6|) into Lagrangian ([2]) 
leads, after straightforward calculations, to the following Euler-Lagrange equations for variables uq(z) and vq(z): 

.duo ( 1 \ 4a 



i- 



dz V 16a 7 4a 



- . , t UqVo, (7) 



■ dv (1 q\ 2 2 

l lz~ = 2 + 2r°-4aTT U °- (8) 



The unperturbed single-color solution to Eqs. J7]), ([5]), u — 0, v (z) — u o exp (— (i/2) (1 + q) z), coincides with 
exact solution Q (with m = 0). Further, to investigate the onset of the parametric instability against the excitation 
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FIG. 1: (a) The critical power for the onset of the parametric instability of the single-color beams with topological charges 
(m — 0) and 1 (m = 2), the beams being unstable at TV > TV C . Chains of rhombuses and dashed and dotted curves show, 
respectively, numerical results and the prediction of the variational approximation produced by Eqs. (|10[) . ()13[) . and (|14|) . 
Examples of periodic oscillations of perturbed solutions at TV > N c are shown for m = 0, q = 0, and TV = 3 in (b), and for 
m — 2, q = 0, and TV = 6.5 in (c). Amplitudes S u o,vo(z) and S u i,v2 are defined by Eqs. (|15[) and (|16|) . respectively. 



of the infinitesimal FF perturbation, we linearize Eq. ([7]) and substitute u Q (z) = u' (z) exp (— (i/4) (1 4- q) z), which 
leads to the following equation for the perturbation's amplitude: 

.du' ( 1 1 q\ , 4a f 



An elementary consideration of Eq. ^ demonstrates that it gives rise to the parametric instability at 
[a + l/(16a) — 1/4 — q/4] 2 < 16a 2 /(4a + l) 2 t> 2 , hence the critical value of total power J5]) is 

1 1 g\ 2 /4a + l x2 



Njr=°)=4x[a+-!----±) . (10) 

V 16a 4 4 J V 4 " / 

A similar stability analysis can be performed for the vortical single-color beam ([4]) with m — 2. In this case, the 
natural ansatz for the FF perturbation is taken in the form of the first angular harmonic, i.e., 

u = u[(z)rexp [-ar 2 +i9- (i/4) (3 + q) z] , (11) 

which, after straightforward calculations, leads to the following evolution equation [cf. Eq. (j9])]: 

du\ /„ 1 3 q\ , 64a 2 . 

^ = { 2a + ^-l-V<~WTW^ M • (12) 

Taking into regard expression ([5]), Eq. (|12j) yields the critical power, 

512a 4 ^ + 8a 4 4 



(13) 



The final analytical prediction for the destabilization thresholds is obtained by the minimization of each expression, 
pop and (|13[) with respect to the variation of free parameter a, for given q. In particular, an explicit result of the 
minimization is that, in both cases of m = and m = 2, 

TV c (m=0 ' 2) ^0 at q > I, 
/V c (m=0 ' 2) « 7T [1 + (m/2)] (1 - g) 2 at 0<1- 9 «1. 

The so obtained critical values of TV C are shown in Fig. [IJa) by dashed lines, which approximate the numerical results 
very accurately. Note that Eq. ([T4"]) explains the vanishing of TV C at q > 1, which is obvious in Fig. Ufa). 

As concerns the realization of the instability in the optical waveguide, it is relevant to note that its length is finite 
(and usually not very large) in a real experiment. For this reason, the observable instability threshold may be shifted 
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to somewhat larger values of N, as a very small instability growth rate would not be able to manifest itself on a 
relatively short propagation distance. 

Above the instability threshold, i.e., at N > Nc m \ simulations of both Eqs. ([T]) and VA equations, (O and (JHJ, 
demonstrate oscillatory behavior of perturbed solutions. An example is displayed in Fig. [IJb), at m — 0, q — 0, and 
N = 3, for variables which, essentially, measure the amplitudes of the zeroth angular harmonic in the FF and SH 
fields: 



SuO,vo( z ) 



{u (x,y) ,v(x,y)}dxdy 



(15) 



It is observed in the figure that the unperturbed state with the zero FF amplitude is periodically recovered. A similar 
example for m = 2, q = 0, and N = 6.5 is displayed in Fig. HJc), for the integral amplitudes of the first and second 
angular harmonics in the FF and SH fields, respectively: 



S,, 



{u(x,y)e ie ,v(x,y)e 2i6 } dxdy 



(16) 



B. Half- vortices: the beams with topological charge 1/2 (m = 1) 



A noteworthy peculiarity of the single-color beams (j4]) is that they may have m = 1 , which formally correspond to 
the topological charge m/2 = 1/2. Of course, this is only possible due to the fact that the FF field is absent in the 
stationary solution; nevertheless, it is shown in what follows below that the half-integer charge of the unperturbed 
solution essentially affects the perturbed evolution of the single-color vortex above the instability threshold. 

To test the stability of these half-vortices within the framework of the VA, a natural ansatz for the FF per- 
turbation is defined as a combination of the zeroth and first angular harmonics: u = u' exp (—a 
u\r exp (— a'r 2 + iO — i'j'z), with 7 and 7' related by the necessary matching condition, 



7 + 7' = 1 + 9/2. 

The VA gives rise to linear coupled equations for perturbation amplitudes u' and u[ [cf. Eqs. © and ([12 

4a 



(17) 
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A straightforward analysis of Eqs. (fT 



[2 (a + a') + 1] 
yields the following resolvability condition: 
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(18) 



(19) 



Substituting relation p7|) in Eq. (fT5|) . one obtains a quadratic equation for 7. The parametric instability sets in at a 



critical value of the total power, N = N, 
to the following result: 



(m=l) 



when the discriminant of the quadratic equation vanishes, which leads 
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(20) 



The expression on the right-hand side should be minimized with respect to the variation of a and a'. The so generated 
critical (dashed) curve vi™ 1-1 ' (q) is displayed in Fig. Ufa) along with the results produced by direct simulations of 
Eqs. ([1]). Like in the situation displayed for m — and m — 2 in Fig. [Ha), the variational prediction is very close to 
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FIG. 2: (a) The critical power for the onset of the parametric instability of the single-color semi-vortex with topological charge 
1/2 (m = 1). Chains of rhombuses and the dashed and dotted curves show, respectively, numerical results and the prediction 
of the variational approximation produced by Eqs. (|20|) and <|21[) . (b) An example of quasi-chaotic oscillations of an unstable 
perturbed solution is shown for m = 1, q — 0, and N = 5.8. (c) The 3D profile of the FF component of the same solution at 
z = 1000. 



its numerical counterpart. Further, also similar to what was done above, an explicit result can be obtained from Eq. 
(j20f in the following form [cf. Eq. (fT4l) ]: 

A c (m=1) =0 at q> 1, 
iV c (m=1) w 2tt (1 - qf at 0<1- 9 «1, ' 

see the dotted line in Fig. HJa). Incidentally, it coincides with Eq. (fT4"|) with m = 2. 

At iV > Nc n ~ 1 \ the perturbed evolution demonstrates persistent quasi-chaotic oscillations of the FF and SH 
amplitudes, on the contrary to perfectly periodic oscillations in the cases of m = and 2, cf. Fig. QJb,c). A typical 
example is displayed in Fig. 2(b) in terms of the integral amplitudes defined as per Eqs. (fT5|) and (fToj) . for q = 
and A = 5.8. In particular, it is observed that the instability generates the zeroth angular harmonic in the SH field, 
represented by amplitude S v q, which was absent in stationary solution that included solely the second angular 
harmonic. Figure [^c) displays, for the same solution, a 3D plot of \u(x,y)\ at z — 1000. The groove structure 
observed in the plot is explained by the fact that \u (x, y) \ — \u' (z) exp(— i~fz — ar 2 ) + u' 1 (z)r ex-p(~ij'z — a'r 2 + i8) | 
becomes small near a certain value of angle 9. The dynamics of the groove is further illustrated in Fig. [3] by a 
set of three contour plots of \u(x,y) |, plotted at z = 1000, 1002 and 1004. It is observed that the groove rotates 
counter-clockwise. This feature resembles the above-mentioned mixed screw-edge dislocation carried by the beam 
with the half-integer vorticity [l3j], although the amplitude does not vanish in the groove, and the inspection of the 
respective phase field does not feature a clear jump by tt, which may be explained by the fact that the quasi-chaotic 
dynamics stirs the phase structure. 



IV. THE STABILITY OF TWO-COLOR BEAMS 

While the single-color beams, with u — 0, are unstable at N > Nc \ in precisely the same region there appear 
two-color modes for m = and m = 2 , built of the FF (u ^ 0) and SH fields. In other words, A^™' is not only the 
stability border for the trapped single-color modes, but also the existence threshold for their two-color counterparts. 
It is relevant to stress that, as seen from Fig. Q]and Eqs. flDij) and (|2"Tj) . the threshold vanishes at q > 1. 

It is easy to construct the two-color solutions for m — 0, in the form of u(x,y,z) = e~ lllz uo(r), v(x,y,z) — 

e~ 2l ^ z v$(r). These solutions exist at A^ > JVc™ -0 ' as completely stable modes, which is natural, as 2D fundamental 
solitons in the system are stable even in the absence of the trapping potential @, Q. 

The transition from the solutions with u — to u ^ 0, with the simultaneous destabilization of the single-color 
state, is a pitchfork bifurcation, which happens with the increase of N, as is shown in Fig. 2] for m ~ 2 (the pitchfork 
must give rise to two mutually symmetric modes at N > N c , which corresponds to the obvious fact that any solution 
with u ^ has its counterpart with the same v and u — > —u). Of course, such an transition does not occur for 
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FIG. 3: Contour plots of the FF component of the same solution which is displayed in Figs. [2jb,c), at z = 1000 (a), z = 1002 
(b), and z = 1004 (c). 
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FIG. 4: Amplitude it ma x of the FF field (a) and propagation constant (b) of the solitary-vortex beams with m = 2 at q = 0, 
of both the single- and two-color types, at N < Nc™ 1 ^ and N > iv| m— 2 ', respectively. The pitchfork bifurcation occurs 
at N = JVi m ~ 2 '. Chains of rhombuses and the dashed curve show, respectively, numerical results and the prediction of the 
variational approximation based on ansatz (122 p . 



the single-color half-vortices with m — 1, as in that case the FF field, u(x,y,z) ^ 0, would carry intrinsic vorticity 
rn/2 = 1/2, which is impossible. 

For the vortex with m = 2, the solution extended past the bifurcation point can be looked for as per the ansatz 



u (x, y, z) — u\qt exp (—ifiz + it 
whose total power is [cf. Eq. ((SJ)] 



or 



! ) , v (x, y, z) = V2or 2 exp (—2i/j,z + 2i8 — /3r 2 ) 



N = 7T 



(2a) _ 2 ul +p- 3 v. 



20 



(22) 



(23) 



with variational parameters uiq, V20 and a, /3. Figure 0] shows the amplitude of u±o and propagation constant \x 
for the vortex beams with m = 2 at q — 0, as predicted by the VA on the basis of this ansatz, and as found from 
numerical solutions. The dashed line in Fig. S^b) represents the constant value, fi = 3/4 + q/2, for the single-color 
vortex with m = 2, which is unstable at N > Ni m ~ 2 \ The VA is very accurate close to the bifurcation point, showing 
a discrepancy which slowly increases at large N . 

A crucially important issue is the stability of the vortex beams at N > Ni m ~ 2 \ as all such states are unstable 
against azimuthal splitting in the free space In the analytical form, the stability can be explored by mean 
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of the nonstationary version of the VA, based on the following ansatz, which adds perturbations containing spatial 
harmonics with numbers —1 and 3 in the FF field, and zeroth and fourth harmonics in the SH component, to the 
stationary solution taken as per Eq. (|22|) (such perturbations induce the splitting instability of the solitary vortex in 
the free space): 



u (x, y, z) — ui(z)r exp (i9 — ar 2 ) + u_i(z)r exp (—id — a_ 1 r 2 ) + u 3 (z)r 3 exp (3i9 — a 3 r ) , 
v (x, y, z) = v 2 (z)r 2 exp (2i0 - (3r 2 ) + v (z) exp (-/3 r- 2 ) + w 4 (z)r 4 exp (4i6> - /3 4 r 2 ) . 



(24) 

Here a and /? are the same as found from the stationary version of the VA based on ansatz (j2"2"]) . Straightforward 
calculations lead to the following linearized evolution equations for the perturbation amplitudes: 
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16/34MiU 3 
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(25) 



w (o) e - i7o ^ ^ = w (o) e -i 7 4^ wit h 

into 



Solutions to Eqs. (|23|) are looked for as u_i 

propagation constants subject to the matching conditions, which ensue from the substitution of expressions 
Eqs. ((251): 

7o = [t- + 7-1)73 = 2/i - 7-1,74 = 3/i - 7_i. 



(26) 



Finally, the eigenvalue problem for 7_i amounts to the following equation, in which 70, 7 3 , and 74 should be substituted 
as per Eq. 



Qa^a^ (8^20) 
(a_i +a 3 +/3) 8 

4a 3 



2a. 
1 

4a 3 
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8a_i 



7-i 



32a 2 . 1 /3 u 2 



(a + a_i + /3o) 4 (4A) + /3c7 + 4? - 870) 



73 



2a 4 /3 4 5 (64ui ) 



(a + a 3 + /3 4 ) 10 (20/? 4 + + 4q - 874) 



(27) 



The critical value of the total power [see Eq. (|23|) ] at the onset of the instability of the trapped vortex ([22j) . TV = 7V C 2, 
is identified as the one at which eigenvalue 7—1, found from Eq. (|27p . becomes complex. Then, it should be minimized 
by varying the set of parameters a_i,a 3 ,/3 and /3 3 , cf. Eqs. (QjJJ), (|T3|) . and (|20|) . Actually, this procedure is too 
cumbersome in its full form, but we have found that a nearly- minimum value of N c2 corresponds to a_i = a 3 = 1.2a 
and /?o = /3 3 = 0.5. The result is reported in Fig. [5](a), which displays N c2 versus q, along with the previously found 
critical value, N c , which is the border between the stable single-color vortices and emerging stable two-color ones. As 
seen from the figure, the two-color vortex is completely unstable at q < q min w —0.5, and has an expanding stability 
area at q > <? m i n . The increase of N c2 with mismatch q is explained by the fact that large q corresponds to the 
cascading limit [H-0j m which the \^ nonlinearity is transformed into the self-focusing cubic interaction, with a 
decreasing effective cubic coefficient, ~ 1/<Z> an d the vortices trapped by the HO potential in the respective weakly 
nonlinear cubic medium have a large stability area (23[-[29j]. The accuracy of the VA prediction for N c2 is essentially 
lower than it was for N c , because the variational ansatze (|22l) and (|24p cannot approximate the complex structure 
of the respective modes accurately enough, and by the above-mentioned fact that the full minimization procedure is 
too cumbersome in this case [in particular, we expect that a more thorough procedure would make the VA-predictcd 
values of N c2 somewhat smaller than those plotted in Fig. [5ja), thus reducing the discrepancy with the numerically 
found stability border]. Nevertheless, the VA predicts the value of q = q m i n almost exactly, and the general shape of 
the stability boundary, N = N c2 (q), is predicted correctly too. 

At N > N c2 , the simulations demonstrate that the instability transforms the stationary two-color vortex into a 
persistent oscillatory one, which keeps the vortical component, mixing it with the zero-vorticity one. To illustrate this 
finding, the peak values (largest over the period of the oscillations) of the integrally defined amplitude of the zeroth 
angular component in the SH field, S v o(z) [see Eq. (|15p], is displayed as a function of N in Fig. [SJb) for q = 0. The 
persistent oscillatory dynamics developed by the unstable two-color vortices at N > N c2 is further illustrated in Fig. 
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(b) 
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FIG. 5: (a) The stability area of the two-color vortices with topological charge 1 (m = 2), in the plane of the mismatch (q) 
and total power (N), is N c < N < N c 2- The numerically found stability borders and their counterparts predicted by the 
variational approximation are displayed by chains of symbols and dashed lines, respectively [in fact, N c (q) is the same border 
as one labeled by m = 2 in Fig. H^a)]. (b) The peak value of the integral amplitude of the zeroth angular harmonic of the SH 
field, Svo(z) [defined as per Eq. (|15[) ]. as found from the simulations of the oscillatory solutions at N > N2c ~ 27.5 for q — 0. 



(a) 



(b) 



24.5- 




23.5- 



800 



FIG. 6: The periodic evolution of the integrally defined amplitudes of the first angular harmonic in the FF field (a), and the 
zeroth harmonic in the SH field (b). 



[S]by the evolution of the integrated amplitudes S u i(z) and S v o(z) [defined as per Eqs. (IT51) and (|15[) ]. produced by 
the direct simulations for q = and N = 39, cf. Fig. [TJc) for unstable single-color vortices. In particular, Fig. [BJa) 
implies that the oscillating mode keeps its vorticity. 

In fact, the oscillating vortices undergo periodic splitting and recoveries. An example of this generic dynamical 
regime as shown in Fig. [7] for the same mode whose evolution is presented in Fig. [5] A similar persistent regime was 
found in the 2D model with the self-attractive cubic nonlinearity, above the threshold of the instability of trapped 
vortices (see details in Ref. |27|). However, in the case of the cubic equation the splitting- recovery scenario is replaced, 
at still larger N, by the onset of the collapse, which does not occur in the x system. 



V. CONCLUSION 



The objective of this work is to demonstrate the possibility of the stabilization of x^ solitary- vortex modes by the 
isotropic trapping potential. The existence and stability of various modes supported by the system is determined by 
the total power, N, mismatch q, and modal topological charge, m/2 (m is the intrinsic vorticity of the SH component). 
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FIG. 7: The regime of periodic splittings and recoveries of an unstable two-color vortex (the same one whose evolution is 
presented in Fig. [6]) is illustrated by a sequence of profiles of \u(x,y,z)\ at z — 90 (a), 140 (b), and 190 (c). 



Using, in parallel, numerical solutions and static and dynamical versions of the VA (variational approximation), we 
have found that, at N < Nc" l \q), all modes are of the single-color type, represented by the SH (second- harmonic) 
component only. In contrast with the usual assumption that such modes are subject to the parametric instability 
against the generation of the FF (fundamental-frequency) field, we have demonstrated that they are stable at N < 
Nc m \q), including the (formal) semi-vortex with topological charge 1/2 (m = 1). However, Nc m \q) = at q > 1, i.e., 
the single-color modes are indeed completely unstable at large q. The modes with m = and m — 2 undergo pitchfork 
bifurcations exactly at N = Nc 7n \ which, destabilizing the single-color states, give rise to stable two-color complexes. 
For m = 0, the emerging states are always stable, while the vortical two-color mode, with m = 2, has an upper 
stability limit, N — N C 2(q). At N > N C 2(q), the unstable vortices feature periodic splittings and recoveries, keeping 
their topological charge. In addition to optics, these results may be realized in atomic-molecular BEC mixtures. The 
semi- vortex does not bifurcate at A = Ai™ -1 '*; instead, it develops persistent quasi-chaotic oscillations, involving 
additional angular harmonics in both the SH and FF fields, and features a rotating "groove" , which resembles the 
mixed screw-edge dislocation induced by the semi-integer vorticity. 

This work suggests possibilities for the analysis in other directions. The stability of trapped vortices and semi- 
vortices with higher values of the topological charge (m/2 = 1.5, 2, 2.5, 3, ...) may be a natural generalization, as well 
as the consideration of trapped modes in the framework of the Type-II (three-wave) x system. On the other hand, 
it may be easy to perform a similar analysis for fundamental and higher-order odd and even (spatially antisymmetric 
and symmetric, respectively) modes in the ID version of the x^ system. In particular, the trapping potential may 
have a chance to stabilize the ID odd modes, which are always unstable in the free space (35[. On the other hand, 
a challenging problem is to find stable solutions for 3D "bullets" in the extended version of system JT]), including 
the temporal variable and group-velocity-dispersion terms. The possibility of the existence of stable "bullets" with 
embedded vorticity is an especially intriguing issue (cf. Ref. (36|, where stable "spinning bullets" were found in the 
free-space model combining the x^ and self-defocusing x' 3 ' nonlinearities). The 3D setting makes it also possible to 
study collisions between stable solitons moving along the trapping potential pipe. 
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